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Abstract. We study the relation between the density profiles of dark matter halos and their 
mass assembly histories, using a statistical sample of halos in a high-resolution N-body simulation 
of the ACDM cosmology. For each halo at z = 0, we identify its merger-history tree, and 
determine concentration parameters c v i r for all progenitors, thus providing a structural merger 
tree for each halo. We fit the mass accretion histories by a universal function with one parameter, 
the formation epoch a c , defined when the log mass accretion rate dlogM/dloga falls below a 
critical value. We find that late forming galaxies tend to be less concentrated, such that c v \ r 
"observed" at any epoch a Q is strongly correlated with a c via c v ; r = c\a /a c . Scatter about this 
relation is mostly due to measurement errors in c v i r and a c , implying that the actual spread in 
c v ; r for halos of a given mass can be mostly attributed to scatter in a c . Because of the direct 
connection between halo concentration and velocity rotation curves, and because of probable 
connections between halo mass assembly history and star formation history, the tight correlation 
between these properties provides an essential new ingredient for galaxy formation modeling. 



1 METHOD 

We investigate the connection between halo density profiles and their mass assembly histories, using 
a structural merger tree constructed from a high-resolution N-body simulation of a flat ACDM model 
with Q m = 0.3, h = 0.7 and o~g = 1.0, whose evolution has been simulated with the ART code g|. The 
trajectories of 256 3 cold dark matter particles are followed within a cubic, periodic box of comoving size 
60/i -1 Mpc from redshift z = 40 to the present. We use distinct halo catalogs at 36 output times spaced 
between z = 7 and 0. NFW density profiles |J, /) NFW (r) = p s /[(r/R B ) (1 + r/R s ) 2 ], are measured for 
each halo with more than 200 particles, corresponding to halos more massive than 2.2 x 10 n h~ 1 M Q . 
For each halo at z = 0, we identify its full merger history, and determine concentration parameters 
c v ; r = R v i r /R s for all progenitors, thus providing a structural merger tree for each of ~ 3000 halos. 



2 RESULTS 

Figure [l] (left) shows the history of mass growth for the major progenitors of several different halos, 
spanning a range of masses and concentration parameters. Massive halos tend show substantial mass 
accumulation up to late times, while the growth curves for less massive halos tend to flatten out earlier. 
By examining a range of full mass assembly histories for our sample of halos, we find that both average 
mass accretion histories and mass accretion histories for individual halos are well characterized by a 
simple function: M (a) = M Q e~ az , a = (l+z)" 1 . Fits to this equation are shown in Figure [l] (left) for 
representative individual halos. The single free parameter a can be related to a characteristic epoch 
for formation, a c , defined as the expansion scale factor a when the logarithmic slope of the accretion 
rate, dlogM/dloga, falls below some specified value, S; the functional form implies a c = a/S. The 
same formation epoch can be defined equivalently for any "observing" epoch z Q of that halo, by 
replacing a by a/a Q , in which case a c = a Q a/S. Thus at any such observing redshift, with scalefactor 




Figure 1: Left: Selected mass accretion trajectories, showing the evolution of the most massive 
progenitor for individual halos in the simulation (thick). Functional fits to the growth curve of each 
halo using Eq. |l]are shown as thin smooth lines. Right: Concentration versus scaled formation epoch 
a c /a Q , for halos at z = (triangles), z = 1 (circles), and z = 2 (squares). The 4 panels correspond to 
different mass ranges. At all masses and redshifts, the concentration parameter c v ; r is well fit by the 
functional form c v ; r = c\a c /a , where c\ ~ 4.1 (represented by the solid line in each panel). 



a = 1/(1 + z ) and mass M Q = M(z Q ), the mass growth is fit by 



M(a) = M Q exp 
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This implies that, for any halo whose mass accretion trajectory follows this functional form, the 
characteristic formation time is the same regardless of the redshift z Q at which the halo is observed. 

We find that the concentration of a halo, defined as c v i r = Rvir/Rei is tightly correlated with the 
characteristic formation epoch as defined above, and that this relation holds at all redshifts when, 
properly scaled by a Q : 

Cvir = cia D /a c , (2) 

where c\ is the typical concentration of halos whose formation time is at the time of measurement, a c = 
a Q . Figure [l] (right) shows that this formula provides a good description of the observed correlation 
between concentration and formation time for halos at all masses and redshifts. As mentioned, the 
typical formation time is a function of mass, but there is significant scatter in a c for a fixed mass. The 
relation defined by Eq. ^ is able to account for the complete mass and redshift dependence of c v i r , and 
for the scatter in c v ir measured for fixed mass halos. This correlation is has important consequences 
for galaxy formation modeling. For further details on this work, see 0]. 
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